Luttinger Liquid at the Edge of a Graphene Vacuum 
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We demonstrate that an undoped two-dimensional carbon plane (graphene) whose bulk is in the 
integer quantum Hall regime supports a non-chiral Luttinger liquid at an armchair edge. This 
behavior arises due to the unusual dispersion of the non-interacting edges states, causing a crossing 
of bands with different valley and spin indices at the edge. We demonstrate that this stabilizes a 
domain wall structure with a spontaneously ordered phase degree of freedom. This coherent domain 
wall supports gapless charged excitations, and has a power law tunneling I — V with a non-integral 
exponent. In proximity to a bulk lead, the edge may undergo a quantum phase transition between 
the Luttinger liquid phase and a metallic state when the edge confinement is sufficiently strong 
relative to the interaction energy scale. 

PACS numbers: 73.43.-f,73.43.Jn, 73.20.-r, 73.22.Gk 



Introduction - Two-dimensional carbon sheets [l|, 
known as graphene, are emerging as one of the most 
exciting new systems supporting the quantum Hall ef- 
fect. This material is different than more standard two- 
dimensional electron gases (2DEGs) because in the ab- 
sence of a magnetic field, the single particle spectrum is 
linear in the vicinity of two inequivalent points in the 
Brillouin zone. The low-energy states near these points 
are described by the Dirac equation j2|i ^.nd in a strong 
magnetic field, the quantum Hall steps that emerge are 
shifted relative to standard 2DEGs [1 H- This effect 
arises because the spectrum of the Dirac equation with a 
magnetic field has doubly degenerate Landau levels (LL) 
for each spin, with one pair at zero energy, half of which 
are filled in the nominally undoped situation. Thus fill- 
ing the Landau states to reach the first quantum Hall 
plateau requires only half as many electrons as needed to 
reach the subsequent plateaus, shifting the step pattern. 




In an undoped standard 2DEG system, there is lit- 
tle interesting electron physics because the filled valence 
states are far below the chemical potential in practical 
situations. By contrast, the partially filled LL at zero 
energy in graphene allow for interesting low-temperature 
physics even in this nominal "vacuum" . When interac- 
tions are included, the half-filled zero energy states rep- 
resent a multicomponent system, which in the absence of 
spin or valley splitting potentials spontaneously polarizes 
due to exchange |5j. In practice, it is spin-polarized in 
the bulk due to the Zeeman splitting. As we show below, 
the vacuum is thus a quantum Hall ferromagnet, with an 
associated low-energy spin-wave. 

In addition to these surprising bulk properties, 
graphene also has an unusual edge structure even in a 
non-interacting picture: the lowest LL has half as many 
edge states as the higher LLs jOj. In this work, we 
demonstrate a remarkable effect when the edge structure 
and the quantum Hall ferromagnetism are both taken 
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FIG. 1: Example of a domain wall configuration. See text. 



into account. Under appropriate circumstances undoped 
graphene forms a a coherent domain wall (DW) between 
the spin polarized state in the bulk and an unpolarized 
region at the edge. The low-energy theory of this DW 
has a U(l) symmetry with a Luttinger liquid Hamilto- 
nian Q- More specifically, the DW may be described by 
a variational wavefunction of the form 



i*>=n 



X<L 



^+sin^e^^Cl,,C,Tx)|Vac>, 

(1) 
create electrons in the two rele- 



where C\\x ^^^ ^-\x 
vant levels, X is the guiding center quantum number with 
allowed values up to the edge at L, |Vac > denotes the 
bulk undoped graphene state (i.e., vacuum) which is par- 
tially polarized since the two spin up lowest LLs are fully 
occupied, and 0{X) and (p are variational parameters. 
An example of 9{X) found by minimizing the effective 
energy functional is illustrated in Fig. 1. The energy of 
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FIG. 2: {Color online) Energy bands for electrons in a 
graphene ribbon with armchair edges in a magnetic field from 
tight binding calculations, ao = lattice constant, unit of en- 
ergy £i = \/27tso/^. ky — X/l"^ with X the guiding center 
coordinate. B=100T, ribbon width is 460A. 



the state is independent of (/), indicating a spontaneously 
broken symmetry in the DW groundstate, with an asso- 
ciated gapless collective mode which may be understood 
as states in which </) has a spatial gradient [a, IE bi|. 
Gradients in (p carry a charge density, and a full 2tt rota- 
tion contains a single electron above the vacuum jlllll2l |. 
Thus this system supports gapless charged excitations. 

This coherent DW may be probed by tunneling into 
it from a bulk metallic lead |jl3|. For a standard 2DEG, 
in the undoped case the system trivially behaves as an 
insulator, and for integer quantum Hall states one finds 
a metallic response [l3|. For the coherent DW, we ex- 
pect a power law tunneling I — V , with exponent deter- 
mined by the pseudospin stiffness and the strength of the 
confinement potential. This is quite different than stan- 
dard quantum Hall edges, where the tunneling exponent 
is thought to be set by the bulk filling factor |lj|. By 
varying the strength of the electron-electron interaction 
(for example, by a screening gate) or the edge confine- 
ment potential, one can vary the tunneling exponent, and 
in in principle may drive the system through a quantum 
phase transition in which the tunneling perturbation be- 
comes relevant in the renormalization group sense. This 
presumably drives the system into a metallic state with 
a linear I —V . 

We now provide details of these results. 

Pseudospin Ferromagnetism in Graphene - The non- 
interacting spectrum of a graphene ribbon in a magnetic 
field with armchair edges is illustrated in Fig. 2 '(^ . The 
eigenstates are characterized by a quantum number X 



which specifies the center of the wavefunction along the 
direction perpendicular to the edge, and spin {a = ±^) 
and valley indices (r = ±1) for a total of four degener- 
ate Landau bands inside the bulk of the system. In the 
undoped system, all the iiegative energy states and two 
of the zero energy states [l^ are filled at zero tempera- 
ture. In what follows, we will ignore the LLs well away 
from zero energy, since these are either completely filled 
or empty and will not affect the physics described here. 
Retaining just the four lowest Landau levels (LLLs) 
near zero energy, apart from constant terms the Hamil- 
tonian may be written as 

H = Y. [-E.<yClxCr.X + A{X)TCl^CraX 
T,a,X 

^ is T. ^'~y^PrA-^)Pr'..'{ci). (2) 

cr,a"',r,r',q 

In this expression, Ez — gfJ-sB is the Zeeman cou- 
pling, A{X) is the energy splitting produced by the edge 
|6| , N^/S is the number of fiux quanta per unit area pass- 
ing through the system, and Vq = ^^2^- is the Coulomb 
interaction. Note that we have assumed an SU(4) sym- 
metric form for the interaction. Although not exact this 
should be an excellent approximation for LLL states pro- 
vided the magnetic length is much larger than the lattice 
constant "2] . Finally, the LLL density operators have the 
form p.,;(q) = j^ Ex e-*«==(2X+,„)c+ ^c'.^x+,„. 

The vacuum state (i.e., undoped groundstate) involves 
filling two of the four LLL. For A = £'2 = 0, any two 
orthogonal states will yield the same energy. In the bulk 
(where A = 0) the Zeeman coupling breaks this symme- 
try, and in the Hartree-Fock approximation the vacuum 
state may be written in the form 



iVac >= H^Ci 



TX "-^-TXl 



<0), 



(3) 



k^ 



where In < O) denotes the state with all the negative 
energy levels filled. Our vacuum is thus ferromagnetic 
The energy per particle of the ground state is —E; 
with S = — \/f 77 the electron self-energy. Note that 
this energy scale is much larger than E^-, and that the 
primary role of E^, is to pick out which two states are 
actually occupied. 

The ferromagnetic vacuum supports four collective 
neutral excitations, which may be approximately gener- 
ated by applying the operators 



'i,TiqJ 






e-i..(2X+.„)C+^C'.,TX+,„ (4) 



to |Vac >. Each of these involves a spin fiip, and two of 
them also include a valley density wave. Because of the 
SU(4) symmetry of the interactions, all four excitations 
are degenerate with energy at long wavelengths iOQ{(i) ~ 



2E^ + Anp sq^i^ , with p^ = l/16V27r in units of e^/eo^, 
and £ = y^c/eB is the magnetic length. This is precisely 
the form one expects for a ferromagnet in the presence 
of a Zeeman coupling. 

It is interesting to notice that due to the spin-charge 
coupling inherent in a LL, doping the system could de- 
grade the spin polarization [l3|. Previous estimates [l3| 
suggest that E^ is too large in the bulk for this to occur. 
However, near the edge where the non-interacting levels 
cross, we now demonstrate that such effects are crucial 
for obtaining a correct description of the lowest energy 
states of the system. 

Coherent Domain Wall at the Edge - Fig. 2 illustrates 
the behavior of the non-interacting levels with armchair 
edges. Because of the Zeeman splitting, a hole-like spin- 
up state necessarily crosses an electron-like spin-down 
state la, [13 • We note that the t = +,— states near 
the edge do not represent states purely in either valley, 
as they are admixed by the boundary condition for an 
armchair edge |a|- Because of its effective spin stiffness, 
the groundstate does not adhere to the lowest energy 
single particle state, but rather generates a DW to make 
the spatial transition. Eq. Q] describes such a DW when 
0{X) passes from to tt. The energy of this state can be 
shown to have the form 12 fl I2II 



E ~ nfp, 



V f— 



X<L 



J2 [AiX] - E,]cose{X), 



X<L 



for an edge at a; = L, provided 0{X) does not vary rapidly 
on the scale of £. This sine-Gordon-like energy functional 
may be straightforwardly minimized numerically |20j . A 
typical example of a DW solution is illustrated in Fig. 1. 
Because the energy of |^ > is independent of 4>, it is 
immediately clear that this represents a broken symme- 
try state which must support a gapless mode. The DW 
is effectively a one-dimensional easy-plane ferromagnet, 
which may be described by an effective action 



S'o= drdy 



d<l)_ 
dy 



-m{y,Tf+- ( ^ 1 4^m(y,r)( ^ 



a7 



(5) 



where r is imaginary time. The field m has its origin in 
the variable 9, and may be understood as the position of 
the center of the DW. The constants F and p may be esti- 
mated from 9{X) and E2 — A{X) using spin-wave theory 
[l0|,|20j. Details of this will be published elsewhere. 

Because 0{X) passes from to tt over the width of 
the domain wall, any configuration of ip{y) that passes 
from to 2tt represents a state in which the sphere of 
allowed directions of an effective spinor constructed from 
the two states (+, t) and (— , J,) is covered precisely once. 
In this situation the state carries an excess or deficit of 
one electron relative to the groundstate |ll[ Il2| | . More- 
over, the twist in cj) may be spread out over the entire 
length of the domain wall, so that charge may be intro- 
duced or removed with arbitarily low energy, in contrast 
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FIG. 3: {Color online) Geometry for tunneling into the co- 
herent domain wall. 



to the gap present for charged excitations in the bulk. 
As we now demonstrate, this has important implications 
the tunneling into the edge of graphene vacuum. 

Tunneling from a Metallic Electrode - The geometry 
for tunneling into the DW is illustrated in Fig. 3. We 
model the lead as a simple Fermi liquid, with action 

ujji,n^kz ,kx 

where /3 is the inverse temperature and p the chemical po- 
tential. The Grassman variable a„fc^ fe^ corresponds to the 
destruction operator for a lead electron in the nth Lan- 
dau level state $„fc t, = i^e*'=«2^+*'''->„(a; - kyi^) 

with Lpn a harmonic oscillator state, and En{kz) — 
[n -\- l/2)a;c + A:^/2m, Wc = eB/mc being the cyclotron 
frequency. The tunnel coupling between the lead and 
reservoir takes the form 



^ttj.n — 



■p^^ dyXnky {y)Kkyk, {^n)i'{y: ^n) + h.C. 
^n nky 



where ■0 is the annihilation operator for an electron in 
the domain wall, which we have taken to be located at 



x = z = 0, and Xnkyiv) 



e'^^yyipn{-kyf). Us- 



ing standard bosonization 7] , one may write the fermion 

operator as ^{y^r) ~ e-'^(^'-)e""^-== '*''"'"''"'. After 
tracing out the lead degrees of freedom, one finds the 
partition function function may be written in the form 

with 



1 






^-^-- 






dyidy; 



Xnky {yi)Xnky (2/2)'0*(2/l, Wn)'0(y2, Wr, 
iUn - p + Enik^) 



Without qualitatively changing the result, one can re- 
strict the LL sum to just the LLL. By noting that 



J2k XX* is strongly peaked at j/i =y2, S may be recast 
in tlie form 

S^-t^ f dTidT2 f dyr{y,Ti)K{n-T2)^j{y,T2). (6) 

Taking the zero temperature limit, one may show J2flj 
K ^ 1/(ti - T2) for large |ti - T2I. 

Our first question is whether the presence of S qualita- 
tively affects the state of the system; i.e., is it a relevant 
operator. A perturbative renormalization group (RG) 
analysis may be applied to S |2(]j , leading to the result 



de 

di 



2)i2 



with the anomalous dimension k = [x + \/x)/2, and 
X = 47r-y//5/r. Estimates of p and F using a spin- wave 
approach ^ yield k « 6.8, 6.0, and 5.3 for B =15T, 
25T, and 45T, respectively. This indicates that under 
usual conditions, S is irrelevant, and the DW remains 
in a Luttinger liquid phase in spite of the coupling to 
the lead. The physical reason for this is that the pin- 
ning energy of the domain wall is small compared to the 
stiffness of the phase angle, because the Zeeman energy 
which ultimately creates the DW is quite small compared 
to the electron-electron energy scale. This suggests that 
enhancing the relative pinning energy can drive the sys- 
tem into a state in which S is relevant, perhaps by ju- 
dicious use of a gating geometry. In this situation the 
coupling to the metallic lead becomes important in the 
low-energy physics, and presumably a current injected 
into the domain wall will behave metallically. 

The irrelevance of S under ambient conditions indi- 
cates that we can treat it perturbatively. The tunnel- 
ing current is then a convolution of the spectral func- 
tions for the lead and the domain wall, separated in fre- 
quency by eV ^ with V the potential difference between 
the two systems [2l|23|. The latter spectral function is the 
Fourier transform of the domain wall correlation function 
< '0(2/,r)-0t(y,O) >~ 1/t'^. The power law dependence 
of this Green's function leads to an anomalous power law 
I — V , a. well-known signature of Luttinger liquids. It 
should be emphasized that this differs considerably from 
edge state tunneling in standard 2DEGs, where the ex- 
ponent is set by the bulk filling factor. The continuously 
varying exponent we find for the edge of the graphene 
vacuum are a consequence of the ferromagnetic nature of 
the undoped state, and its unusual crossing of particle- 
and hole-like edge states. 

In summary, we have studied the armchair edge of 
graphene in its undoped state. We demonstrated the 
ferromagnetic nature of the system due to interactions, 
and showed how this leads to a coherent domain wall 



due to a crossing of the non-interacting edge states. The 
domain wall can be probed by coupling it to a normal 
lead, and we found this supports a power-law tunneling 
current and the possibility of a quantum phase transition 
into a metallic state. 
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